Matrica mase

m0
0md

Matrica prigusenja

c + cd-cd
-cdcd

Matrica krutosti

k + kd-kd
-kdkd

Vektor sila

F0*sin(p*t)
0

Pretpostavljeno rešenje u

A*sin(p*t) + B*cos(p*t)

Pretpostavljeno rešenje ud

C*sin(p*t) + D*cos(p*t)

Grupisanje uz sin u prvoj jed. kretanja

A*k + A*kd - A*m*p**2 - B*c*p - B*cd*p - C*kd + D*cd*p - F0

Grupisanje uz cos u prvoj jed. kretanja

A*c*p + A*cd*p + B*k + B*kd - B*m*p**2 - C*cd*p - D*kd

Grupisanje uz sin u drugoj jed. kretanja

-A*kd + B*cd*p + C*kd - C*md*p**2 - D*cd*p

Grupisanje uz cos u drugoj jed. kretanja

-A*cd*p - B*kd + C*cd*p + D*kd - D*md*p**2

Sistem jednacina

A*k + A*kd - A*m*p**2 - B*c*p - B*cd*p - C*kd + D*cd*p - F0
A*c*p + A*cd*p + B*k + B*kd - B*m*p**2 - C*cd*p - D*kd
-A*kd + B*cd*p + C*kd - C*md*p**2 - D*cd*p
-A*cd*p - B*kd + C*cd*p + D*kd - D*md*p**2

Matrica koeficijenata Ms

k + kd - m*p**2-c*p - cd*p-kdcd*p
c*p + cd*pk + kd - m*p**2-cd*p-kd
-kdcd*pkd - md*p**2-cd*p
-cd*p-kdcd*pkd - md*p**2

Vektor slobodnih clanova bs

F0
0
0
0

Determinanta sistema Delta

k + kd - m*p**2-c*p - cd*p-kdcd*p
c*p + cd*pk + kd - m*p**2-cd*p-kd
-kdcd*pkd - md*p**2-cd*p
-cd*p-kdcd*pkd - md*p**2

Determinanta za A

F0-c*p - cd*p-kdcd*p
0k + kd - m*p**2-cd*p-kd
0cd*pkd - md*p**2-cd*p
0-kdcd*pkd - md*p**2

Determinanta za B

k + kd - m*p**2F0-kdcd*p
c*p + cd*p0-cd*p-kd
-kd0kd - md*p**2-cd*p
-cd*p0cd*pkd - md*p**2

Determinanta za C

k + kd - m*p**2-c*p - cd*pF0cd*p
c*p + cd*pk + kd - m*p**20-kd
-kdcd*p0-cd*p
-cd*p-kd0kd - md*p**2

Determinanta za D

k + kd - m*p**2-c*p - cd*p-kdF0
c*p + cd*pk + kd - m*p**2-cd*p0
-kdcd*pkd - md*p**20
-cd*p-kdcd*p0

Rešenja [A=DeltaA/Delta, B=DeltaB/Delta, C=DeltaC/Delta, D=DeltaD/Delta]

F0*(cd*p*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)) + (kd - md*p**2)*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)))/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)F0*(cd*p*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)) - (kd - md*p**2)*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)))/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)F0*(cd*p*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)) + kd*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)))/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)F0*(cd*p*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)) - kd*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)))/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)

Amplituda mase m: u =

sqrt(F0**2*(cd*p*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)) + (kd - md*p**2)*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)))**2/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)**2 + F0**2*(cd*p*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)) - (kd - md*p**2)*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)))**2/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)**2)

Amplituda mase md: ud =

sqrt(F0**2*(cd*p*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)) + kd*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)))**2/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)**2 + F0**2*(cd*p*(cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2)) - kd*(-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2)))**2/((-2*cd*kd*p + cd*p*(k + kd - m*p**2) + p*(c + cd)*(kd - md*p**2))**2 + (cd**2*p**2 - cd*p**2*(c + cd) - kd**2 + (kd - md*p**2)*(k + kd - m*p**2))**2)**2)